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1. Introduction 

Quantities related to nucleon structure are among the most challenging that can be computed 
within lattice QCD. During the past few years, there has been significant progress in these cal- 
culations: quark masses used nowadays are close to the physical ones, lattice spacings are small 
enough to allow for a controlled extrapolation to the continuum limit and volumes are sufficiently 
large to suppress finite volume effects. But still even for some simple observables like the nucleon 
axial charge gA or the average momentum of the unpolarized isovector parton distribution (x) u -d 
the values obtained from lattice calculations differ from their corresponding experimental values. 
In Fig. 1 we show the relative deviation of lattice results for < x > u -d from the value obtained in 
Ref. [1] (ABMK). We also show the deviation of results from different phenomenological extrac- 
tions from deep inelastic scattering experiments as compared to that of ABMK. The deviations of 
the phenomenologically determined values for < x > u -d between the different groups might be in- 
terpreted as a systematic error related to different data sets and different analysis/fitting strategies. 
However, as can be seen, the spread in the phenomenological number is much less than the relative 
deviation between lattice data and phenomenology. One possible reason for the discrepancy can be 
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Figure 1: With the filled symbols we show the deviation of the values of < x > u _ ( / obtained from lattice 
calculations using twisted mass fermions (TMF) [2, 3] (ETMC), a hybrid action of domain wall fermions 
(DWF) on a stagerred sea [4] (LHPC), Clover fermions [5, 6] (QCDSF) and DWF [7] (RBC/UKQCD). The 
open symbols show the deviation of the phenomenological results [8] (JR) and [9] (MSTW) as compared to 
ABMK [1]. 

systematic effects that enter in the lattice calculations. Therefore the latter have to be thoroughly 
studied, which is a current important activity of many lattice collaborations. In these proceedings 
we focus on one potential systematic effect, namely excited state contributions. Furthermore we 
restrict ourselves to the two observables mentioned before, gA and (x) u _ d , that are extracted from 
the matrix element at zero momentum transfer. These quantities are calculated in lattice QCD by 
studying the asymptotic behavior of suitable ratios of 3-point and 2-point correlation functions for 
large Euclidean time separations. 

Several studies of gA and {x) u _ d at different values of the lattice spacing and various vol- 
umes [2, 7, 4, 5] indicate that for the pion masses considered, discretization and finite volume 
effects are likely not sufficient to explain the disagreement of the lattice calculations with the cur- 
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rent experimental value of gA and the value of (x) u _ d from global analyses. Moreover, recent 
preliminary results [5] shows that the discrepancy persists even for pion masses almost as small as 
the physical one. 

In order to clarify whether and to what extent excited state contributions affect the lattice 
calculation of matrix elements, we have performed a dedicated high-precision calculation of the 
3-point functions for gA and (x) u _ d . To be more precise, we have used roughly 7500 measurements 
for gA and about 23000 for {x) u _ d , which enabled us to calculate the correlation functions needed 
for gA and (x) u _ d to sufficient accuracy even for large Euclidean time separations. This allowed us 
to study possible excited state contributions. 

2. Lattice techniques and details 

In this work, we have employed twisted mass Wilson fermions [10] at maximal twist. This 
lattice discretization has the advantage of an automatic 0(a)-improvement, thus removing the ne- 
cessity of additional operator specific improvements. We use gauge field configurations generated 
by the European Twisted Mass Collaboration (ETMC) with Nj ■ = 2 + 1 + 1 flavors and refer to 
Ref. [1 1] for further details of our lattice formulation. 

To understand how the excited state contributions emerge from the calculation of nucleon 
matrix elements the main observation needed here is that to leading order we find two additional 
time dependent contributions for a matrix element (N\<ff\N), namely 

(Af|^|A0 latdce = (N\0\N)+Aexp(-AMt')+Bexp[-AM(t-t')] . (2.1) 

Here AM is the mass gap between the nucleon ground state and the first excited state and A and 
B are constants depending on the particular choice of the lattice nucleon creation and annihilation 
operator, respectively. For details we refer the reader to Ref. [12]. For our precision calculation we 
restrict ourselves to only one ensemble with a pion mass of m % k, 380 MeV and a lattice spacing 
of a « 0.078 fm. This pion mass is chosen sufficiently small in order to be relatively close to 
the physical pion mass on the one hand but still large enough to ensure that finite size effects 
can be safely neglected. Also, for such a pion mass, the calculations of the propagators do not 
require large computer resources, thus enabling us to collect the statistics necessary for the high 
precision we are aiming at. We expect discretization effects to be small due to the automatic 
f?(o , )-improvement, which is confirmed by previous calculations of nucleon matrix elements at 
three different lattice spacings smaller than 0.1 fm [2, 13, 14] using Nj = 2 maximally-twisted- 
mass fermions. Let us point out in particular, that even though this analysis was performed using 
twisted-mass fermions, the important properties of excited state contributions are expected to be 
universal and independent of the particular lattice discretization used. Therefore we expect our 
conclusions to hold qualitatively also for different lattice actions. 

For the precision calculation of the 3-point functions performed in this work we employed a 
method that is different from the standard sequential method, where the sink time slice is fixed. We 
fixed the time slice of the operator insertion as well as the operator to obtain a result for all sink 
time slices. This enables us to study the source-sink separation dependence of a particular matrix 
element. Note that in general the standard calculation is more desirable since it does not depend 
on the operator itself, whereas in the open sink method we have to restrict to a single operator. The 
interested reader is referred to Ref. [12] for a more detailed discussion. 
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As stated in the introduction, we concentrate on two relatively simple but nonetheless phe- 
nomenologically relevant quantities. The first is the nucleon axial charge, gA, playing an important 
role in the beta decay of the neutron and appearing as a low energy constant in effective chiral La- 
grangians. It has been experimentally measured precisely and it is also straightforward to calculate 
in lattice QCD. However, the values obtained from various lattice calculations are typically 5% to 
10% below the experimental result [15] while having themselves a statistical accuracy of the order 
of 1%. The second observable is the lowest non-trivial moment of the unpolarized parton distribu- 
tion function in isovector flavor combination, (x) u _ d . It is determined phenomenologically from a 
global analysis of deep inelastic scattering data, and the discrepancy between the phenomenolog- 
ical and lattice values is even larger, about 50% to 60%, as can be seen in Fig. 1. We would also 
like to stress that non-perturbative renormalization is employed for the bare matrix elements. The 
corresponding renormalization factors are calculated in the RI'MOM scheme and are matched to 
the MS scheme at a scale of (2 GeV) 2 . For more details we refer to Refs. [16, 3]. The values of the 
renormalization constants relevant here are Za = 0.11 A for the renormalization of the bare gA and 
Zm = 0.998 for the renormalization of (x) u _ d [17]. 

3. Results 

We performed an analysis of gA on a single Nf = 2 + 1 + 1 ensemble characterized in Sec. 2 
using the open sink method. The time slice of the operator insertion was fixed to t' = 9a, chosen 
to safely suppress excited state contributions from the source, as can be verified from the 2-point 
function. We have used Gaussian smearing of the quark fields, including APE-smeared gauge links, 
in order to improve the overlap with the nucleon ground state. The result of the analysis using the 
open sink method is shown in Fig. 2 (left). In this figure we also compare the values obtained 
from the open sink analysis to the value of a standard analysis with a fixed source-sink separation 
of 12c? on the same ensemble. The value of gA does not show any dependence on the source-sink 
separation t within statistical accuracy and thus demonstrates the absence of significant contribution 
from excited states. It is worth mentioning that, in order to reach a comparable statistical accuracy 
as the one obtained when using the fixed sink method with t = Via with 500 measurements, we had 
to perform roughly 7500 measurements when we take e.g. t = 18a. We have run the same analysis 
for (x) u _ d . For the open sink method, we have chosen the operator insertion time to be t' = 11a. 
We expect that for this choice excited state effects from the source are sufficiently suppressed. 
With a statistics of 23,000 measurements for {x) u _ d at a source-sink separation of t = 18a, we 
could equal the precision of the fixed sink method that was done with a source-sink separation of 
t = \2a using 1300 measurements. In the right panel of Fig. 2 we plot (x) u _ d as a function of the 
source-sink separation t. We also indicate the value obtained from the fixed sink method analysis 
as well as the experimental result from a global analysis [1]. We observe that the values for {x) u _ d 
reach a plateau at larger values of the source-sink separation than what we have used in the fixed 
sink method and that the plateau value is shifted. Still, despite the fact that the results for larger 
values of t decrease they clearly do not reach the phenomenological value. In order to estimate the 
residual dependence on t, we determined the value of (x) u _ d for an infinite source-sink separation 
by fitting the expected exponential behavior, 

(x) u _ d +Aexp [-AM (t-t')] , 
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Figure 2: Left panel: Results for gA for a range of source-sink separations obtained from the open sink 
analysis on one Nf = 2 + 1 + 1 ensemble. The light gray band indicates the result obtained from the fixed 
sink method using a source-sink separation of 12a and the dark gray band shows the experimental value. 
Right panel: (x) u _ d for a range of source-sink separations obtained by means of the open sink method. 
The value (including errors) obtained from the fixed sink method using a source-sink separation of 12a is 
indicated by the light gray band. The phenomenologically extracted value is shown with the dark gray band. 
The blue solid line corresponds to a fit described in the text. 



to the lattice results with a fixed t' = 11a. The result of this fit is (x) u _ d = 0.22(1). This value 
is 12% lower than the result of {x) u _ d = 0.250(6), obtained using t = 12a in the fixed sink 
method. We estimated the error of the fit by varying the fit range and by comparing the results 
obtained by using a fixed parameter AM as well as by including AM in the fit. We have also 
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Figure 3: Ratio of 3- to 2-point functions for gA (left) and (x) u _ c / (right) versus the time separation of the 
operator insertion from the source. The sink-source time separation — 12a. Nf — 2 gauge configurations 
are used with pion mass ~ 300 MeV and a = 0.089 fm. 

used a generalized eigenvalue (GEV) approach as first suggested in [18] and further developed 
and refined in Refs. [19, 20, 21]. We considered two variational basis: in the first, two inter- 
polating fields are considered, namely the standard one Jn(x) = e abc {d Ta {x)Cysu b {x))u c {x) and 
J' N (x) = e abc {d Ta {x)Cu b {x))y^u c {x). The latter is known to have small overlap with the nucleon 
state but a large one with the Roper. In Fig. 3 we compare the ratio of the three-point to the 
two-point function arising when using the GEV approach to the one using just the standard Jn{x) 
interpolating field as we vary the time separation of the operator insertion form the source. As 
can be seen, for both gA and {x) u -d the ratios are consistent and lead to the same value for these 
observables. In the second, three different levels of smearing are employed to calculated a 3 x 3 
correlation matrix. In Fig. 4 we show the effective plateau matrix elements for fixed ? op = ? s ink/2 
and various to. From the preliminary results as shown in Fig. 4, we conclude that within the statis- 
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tical errors we can achieve with our nucleon 3 -point functions, the values obtained from the GEV 
are fully consistent with our standard method. 



6.5 
6 

5.5 
5 

4.5 
4 

3.5 



-1 


1 










1 


standard 




H 


GEVP, 3x3 — " 



-0.4 



-0.2 







0.2 



0.8 



1 



1.2 



0.4 0.6 

< x > u - d 

Figure 4: Comparison between the value obtained from a standard calculation of (x) u _ d using a fixed source- 
sink separation of 12a and preliminary results from a GEVP analysis using the same source-sink separation 
for a range of ?q. 



4. Summary and Conclusions 

We have performed precision calculations of gA and (x) u _ d for a single ensemble of gauge 
field configurations corresponding to a pion mass of about 380 MeV with Nf = 2 + 1 + 1 dynamical 
fermions employing a non-perturbative renormalization. We have studied the dependence of these 
quantities on the source-sink separation in order to assess the influence of excited states on the 
current lattice results for gA and (x) u _ d . This is particularly important given that excited states 
may play a role in explaining the presently observed discrepancy between lattice computations and 
phenomenological evaluations of several important observables related to nucleon structure. 

We find that for the here considered pion mass of about 380 MeV and lattice spacing of a « 
0.078 fm, the contamination of excited states is negligible for gA, but for (x) u _ d we observe an 
effect that is of the order of 10% compared to our previous calculations, in which the source-sink 
separation has been fixed to approximately 1 fm. This is an effect larger than the finite volume and 
lattice spacing effects found at this value of the pion mass, volume and lattice spacing. Moreover, 
this study shows that excited state contributions are operator dependent and should be investigated 
separately for each operator. 

Motivated by a study where variational methods have been employed in the calculation of 
the B*Bn coupling [21], we applied the generalized eigenvalue method [18, 19] to calculate these 
quantities. The results of this analysis were consistent with the standard method. This indicates 
that in order to extract excited state effects, for which the GEVP method is expected to be very 
efficient, a significantly higher accuracy of the lattice data are needed. 
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